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Abstract. In this paper, we study transmission line resonances in high Tc dc 
SQUIDs. These resonances are exhibited in the characteristics of SQUIDs which are 



T3 

q ' fabricated on substrates with high dielectric constant, such as strontium titanate. The 

O . power balance equation is analytically derived both for symmetric and asymmetric 

SQUIDs. Using this, we investigate SQUID current - voltage I(V), voltage - flux U($) 
and voltage modulation AV characteristics. 
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1 1. Introduction 

Superconducting quantum interference devices (SQUIDs) are the most sensitive 
magnetic flux sensors known today. Studying and improving SQUID properties has 
been a major research field in the last few decades [1, 2}. Most high Tc SQUIDs are 
fabricated on strontium titanate (STO) substrates. It has been seen that large dielectric 
constant of STO strongly affects the performance of high Tc dc SQUIDs by causing 
transmission line resonances to appear in their current- voltage characteristics [3, 4, 5]. 
This was pointed out for instance, by Lee et al who experimentally showed that their 
hairpin shaped SQUID loop could be treated as a quarter wave microwave resonator 
in order to predict the voltage at which the resonance appeared in their measurements 
[6]. Enpuku et al used a theoretical approach to predict the effect of the resonances 
based upon a transformation of the distributed resonator structure of the SQUID loop 
into a series of lumped element resonators [3]. This enabled them to calculate the 
current- voltage characteristics of the SQUID and later, the voltage-flux curves for other 
devices. However, it was necessary for them to solve the resulting differential equations 
numerically In this paper, we develop an analytic approach to investigate the effects 
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of transmission line resonances on dc SQUID resonance characteristics. We are able 
to obtain good agreement between experiment and theory for devices in the literature 
[6]. We find that the resonance positions are mainly controlled by the product of the 
junction critical current 1$ and junction normal state resistance Rs(Io x Rs) as wen as 
the dielectric constant of the substrate er. We also investigate the effects of introducing 
asymmetry in junction parameters. We find that the resonance positions are not very 
sensitive to the asymmetry. The actual magnitude of the current-voltage curve and AV 
curves are however sensitive to the degree of asymmetry introduced. We ignore the 
effects of thermal noise in performing our calculations. Noise analysis can be performed 
perturbatively as in [7]. In such calculations, the lowest order result is the noise free 
case where our analysis will be useful. 

We begin with a brief outline of relevant previous work in the study of the effect 
of transmission line resonances on SQUID characteristics in section 2. In section 3, we 
derive the power balance relation in terms of SQUID parameters, which involves deriving 
an expression for the circulating current through the SQUID inductance. Section 4 
discusses current-voltage /(V) curves for symmetric SQUIDs whereas in section 5, effects 
of introducing asymmetry on I(V), voltage-flux V(<&) curves and SQUID modulation- 
current AV(J) curves are discussed. The paper ends with a discussion of results obtained 
in section 6. 
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Figure 1. (a) Geometry of the SQUID washer. Here I is the length of the slit, w is 
the width of the electrode and s is the width of the slit, (b) Equivalent circuit of the 
washer when its parasitic capacitance distributing along the slit of the washer is taken 
into account. The circuit consists of the SQUID coupled to a series of LCR circuits. 
Rd is a damping resistance in parallel to the SQUID. 
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2. SQUID circuit equations 

The geometry of the SQUID washer that we have used is shown in Fig. (la). This is a 
geometry commonly used to manufacture SQUIDs [6, 8, 9] and is the geometry Enpuku 
et al used to numerically investigate the effects of large dielectric constant of strontium 
titanate (STO) on the characteristics of high T c dc SQUIDs [3]. The slit of the SQUID 
washer makes up the SQUID inductance where /, s, w and d denote the slit length, slit 
width, electrode width and thickness of electrode respectively. For this geometry, the 
inductance per unit length of the slit L and parasitic capacitance per unit length C are 
given by [3, 10]: 

L = L M + L K (1) 

where Lm is the magnetic inductance per unit length and Lk is the kinetic inductance 
per unit length of the SQUID slit given by [11]: 



J M 



K(k') 

2/xqA 2 w 4ws w Aw(2w + s)- 



(2) 



b H-^——t] + 7T— H " , V ]] (3) 



dwk' 2 K 2 (k') s d(w + s) 2w + s d(w + s 



and 



where \x§ is the permeability of free space, K(k) is the complete elliptic integral of the 
first kind [12] with a modulus k = — , k' — (1 — k 2 ) 1 ^ 2 , A is the penetration depth of 

the film, 6r is the dielectric constant of the STO substrate and c is the velocity of light 
in vacuum. The SQUID slit behaves as a transmission line with distributed inductance 
L and distributed capacitance C. The impedance Z^b of the slit seen from terminals A 
and B is given by [3]: 



J AB 



iZ tan(fi/v / LC') + iVtL PR (5) 

where Z = \Jl/C is the characteristic impedance of the slit transmission line, f2 is the 
angular frequency of measurement and Lpr is the junction parasitic inductance. The 
first term arises because the hairpin shaped slit can be treated as a shorted transmission 
line of length I. 
Now, using the formula [3] 



(it x . 4^z< ^ — ^ X 

T ) = "S »»-(2b-1)' 1 (6) 



eqn.(5) can be expanded as (in the lossless case with L PR = 0): 

oo -y 

Zab = 2 MC n + ijim7 n ' (7) 

with C n = CI/ 2 and L n = 8Ll/7r 2 (2n — l) 2 . This transformation to an equivalent cir- 
cuit has the advantage that it allows us to consider loss in the transmission line. In the 
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lossy case, the rf-loss R n is added to eqn.(7) leading to: 



ZAB 5 %nC n + l/iSlL n + 1/Rn ' (8) 

where R n = Q^J L n /C n and Q is a quality factor^. Here / is the SQUID slit length. 
This expression allows us to express the impedance Z^b by the series of L-C-R resonant 
circuits as shown in fig. (lb). 

We will now show derivations of various quantities in the lossless limit and then 
quote the results obtained for the lossy case which can be derived analogously. The 
circuit equations can be easily derived by the application of Kirchhoff's laws. In fig. (lb), 
current entering point C should equal current leaving point C. We are assuming the 
most general case in which the SQUID is made up of junctions which are asymmetric 
[13, 14]. Let the average junction critical current be Jo, the average junction normal 
state resistance be R$ and the average junction capacitance be Cj. Let the asymmetry 
parameter in Jo be k, that in Rs be p and that in Cj be x- So specifically, let 
us split [13] 2J = (1 + k)I q + (1 - k)I ,2/R s = (l+p)/R s + (l-p)/R s and 
2Cj = (1 + x)Cj + (1 — x)Cj. Then, this gives for the first junction: 

(1 + k)I sin(0i) + — + (1 + X)Cj— - — - - = Rd , (9) 

and writing an analogous equation for point D gives for junction 2: 

(1 - k)I sm{6 2 ) + — + (1 - X)Cj— - — + - = ^ , (10) 

Here J is the circulating current through the SQUID inductance, I B is the SQUID bias 
current, V\ and V 2 are voltages across junctions 1 and 2 , 9 2 and Q\ are the phases of 
junctions 2 and 1 and Rd is a damping resistance in parallel to the SQUID inductance. 
The circulating current in Fourier space is given by§: 

V 2 -Vi + %5{tt) , 
J = 2— — — — , (11) 

^AB 

where 5 is the Dirac delta function and $ is the total flux applied to the SQUID loop. 

We will follow [3] and use normalized circuit equations. We normalize currents 
by Jo, voltage by IoRs and time t by $o/2nIoRs- The ac Josephson relation gives 
v = dO/dr = V/IqRs , where v is the normalised voltage and r is the normalised time. 
Here we have neglected random noise currents. The normalized circuit equations are: 

(1 + X )Pch = \(*b + 3) - (1 + p)0i - (1 + «) sin(^) - 7^ - 6 2 ) , (12) 

(1 - x)Pch = \(iB- 3) - (1 - PA - (1 - «) sm(^) + 7(^1 " ^2) , (13) 

\ It is possible to consider different quality factors for each n but we will not do so here. 
§ The expression for circulating current can be written explicitly in the frequency domain in terms of 
Vi and V2 which are obtained by solving eqns. (9) and (10). The equivalent expression for j in the 
time domain (eqn.(B.22)) can only be explicitly obtained after solving the equations and taking into 
account the properties of the SQUID loop. This is presented in Appendix B. 
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where (3 C — 2vrJ Cj-R|/$ is the SQUID McCumber parameter with $ being the flux 
quantum and 7 = R$/Rd- In principle, these equations can be solved numerically 
using the dielectric constant 6r as a fit parameter having a finite number of resonant 
circuits. Enpuku et al present numerical solutions of eqns.(12) and (13) in their paper. 
As shown in the appendix A, they however consider a finite number 'n' of resonant 
circuits and the expression they use for the circulating current J becomes indeterminate 
in the continuum limit. The analytic handle will enable closed form expressions for 
Ib(V), which can subsequently be used to simulate graphs quickly and conveniently. 
Normalizing eqn.(ll) and the expressions for impedance in the lossless limit, eqn.(7) and 
lossy limit, eqn.(8), the expression for the normalized circulating current j in Fourier 
space is given by: 

] = (Q 2 - ~Q X + ^H0)A(u;) (14) 

where u = Q&o/2nI Q Rs and is dimensionless. Here is the externally applied flux 
normalized to <3> and (3 = 2Z// / < ^o is the SQUID inductance parameter. In the lossless 
case with Lpp = 

A (uj) = 2 $0 (15) 
LI 27t/ tan^Tr/o^sVZC/^o) 
and in the lossy case, performing the sum in eqn.(8): 

A(u) = 2i-[tt(C - i£) + *(C + i0\ , (16) 

where ( = 1/2 + itulVW '/2ttQ, £ = iluj^W(l - 4Q 2 )/2ttQ, v = -Q%/2nI VLl, 

/j, = Lnuj^JlC(l — 4Q 2 ) and ty(x) = d x log r(a;) is the di-gamma function with T(x) 
being the standard gamma function [12]. If Lpp 7^ then A(uS) should be replaced by 
A(u)/(1 + nI A(uj)Lp R /&o) in the calculations. 

3. Power balance equation 

The calculation details for the SQUID power balance equation have been included in 
Appendix B. After some straightforward but tedious algebra, the power balance relation 
so derived is as follows: 

i b = 2v dc + ^ C -, 2 ((l + K 2 ){dd + rf] + 2n[df + dr] 



2\ds — rr\ 2 

+ (1 - k 2 ) {[dd- rf] cos(2tt0) +i[dr- dr] sin(27T0)^ | 

+ 2^2-rf^ 1 + 2l ~ 4^ } { (1 + + ^ " + Sf] 

- (1 - k 2 ) {[I3 2 ss - rf] cos(2tt0) - i(3[sf - sr] sin(27T0)^ | . (17) 

Here d = d(vd c ), s = s(vd c ), r = r(vd c ), A = A{vd c ) and x denotes the complex conjugate 

of x. The phase difference between the two junctions have been set to 2n(f). Here, 

2v dc (l + viy/*-v dc (18) 
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s(u) = {oo 2 (3 c + iu) , d(u) = A(u) - {uj 2 (3 c + + 2iury) , (19) 

r(u) = ipuo + xPcuj 2 , (20) 

Note that the power balance relation has the correct symmetries. For instance 
interchanging the two junctions leads to k — > —k, p — > — p, x ~ * ~X an d V'i ^ 4 , 2- The 
expression is symmetric under this transformation. Another interesting observation is 
the fact that the p and x dependence are both packaged in the function r(cu). 



4. Symmetric SQUIDs 

For symmetric SQUIDs, i.e. k = p = x — and in the lossless case where A = A, we 
have 

h = 2vdc + ^Tjtj 1 + cos2vr0| + ^(1 + 2 7 ){l - cos2tt0| , 

(21) 

where s, d are given by eqn.(19) and A is given by eqn.(15). In reference [18], En- 
puku et al derive the SQUID power balance equation in the absence of transmission 
line resonances. By allowing the slit length to tend to zero, hence treating the SQUID 
inductance as a lumped inductance instead of a distributed parameter, we retrieve En- 
puku's expression for the power balance relation: 

. o2 I 2 (l + cos(27r0)) J 2 (l-cos(27r0))(l + 2 7 )^i 

BdC ~ dc+ 2(l + ®vtt + 2([f-^ c F + ^ c (l + 2 7 ) 2 )^ j 

Here we have neglected the capacitance in parallel to the shunt resistance, which they 
had included in their circuit as our analysis is based on a circuit, which does not involve 
such a capacitance term. Lee et al [6] had seen a current step in their /(V) curve 
which they attributed to the transmission line resonances due to the high dielectric 
constant of the substrate. We have used their SQUID parameter values (/ = 55pm, 
SQUID inductance, L SQ (L x I) = 55pH,I = A.75pA,R s = 13.8ft, SQUID parasitic 
inductance L PR = 13pH and j3 c — 0) to simulate ib(vd c ) curves from eqn.17 in the 
lossless case. In accordance with Lee et al, €r = 1930. They had seen the resonant 
current step at 80pV. Fig. (2a) shows the comparison between the curves we get 
from eqn.(17) and eqn.(22) respectively. Fig. (2b) shows the comparison between the 
differential conductance dib/dvd c {vdc) curves obtained from differentiating eqn.(17) and 
eqn.(22) respectively. Enpuku et al's formula does not lead to any resonances which 
clearly arise from the consideration of the SQUID loop behaving as a transmission line. 
The minima in the ib{vd c ) curve occur for A — > oo whereas the maxima occur for A — > 0. 

Eqn.(17) is the SQUID power balance equation when transmission line resonances 
due to the distribution of the parasitic capacitance due to high dielectric constant of the 
STO substrate along the slit of the SQUID washer are taken into account. The reso- 
nances only affect the terms of this equation which contains A(oo). This is because these 
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(a) (b) 

Figure 2. (a) Comparison between the ib{vdc) curve from eqn.(17) which takes 
transmission line resonances into account (solid black line) and the ib{vdc) curve from 
the equation derived by Enpuku et al, eqn.(22) which considers the SQUID as a lumped 
element (dotted line) for = 1/2 using Lee et al's parameters. The step in the solid 
black line at around Vd c — 1-2 represents a resonance position. Such features are clearly 
absent in the dotted line, (b) Comparison between the corresponding dib/ dvdc{vdc) 
curves. The solid black line is derived from eqn.(17) and the dotted line from eqn.(22). 



terms are associated with D = 9\ — 9 2 and hence with the circulating current. Thus, 
where cos27T0 = 1, ie for § ext = n$o the resonances have no effect on the current- 
voltage characteristics, whereas for cos27T0 = —l,$ ext = \{2m + l)<E»o, the resonances 
have maximum effect. 
The maxima occur at: 



Vdc 



m$ 



2H Rs VLC 
The minima occur at: 



(23) 



i> + l)* 1 

4U Rs Vic y ' 



If we approximate VLC by y ^fgr^ , where c is the speed of light in free space, 
then eqn.(24) is the one used by Lee et al to predict the voltage at which the resonance 
occurred for their device. This approximation is only true if the kinetic inductance of 
the SQUID loop is negligible. 



5. Introducing asymmetry 

5.1. Current-voltage characteristics 

In this section we wish to consider the effects of introducing asymmetry to the SQUID. 
We have introduced current and resistance asymmetry systematically and seen the effect 
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on the I(V) curve both for the = and = 1/2 cases. Fig. (3a) shows the comparison 
between Lee's experimental I(V) and our theoretical simulations for the = case 
and fig. (3b) shows the same for the = 1/2 case. The comparison between experiment 




(a) (b) 

Figure 3. (a) Comparison between Lee et al's experimental data read out from 
[6] (bold black dots) and theoretical simulations from cqn.(17) for the = case. 
Continuous line is the simulated curve for no asymmetry considered, dotted line is 
for the case when current asymmetry is considered and k — 0.5, the dashed line is 
for the case when normal state resistance asymmetry is considered and p — 0.2 and 
the dot-dashed line is for the case when both current and resistance asymmetry are 
considered and k — 0.5 and p = 0.2. (b) The 0=1/2 case. 

and theory in the = case is not very good at large voltages. In the 0=1/2 case, 
introducing current asymmetry shifts the origin to a point higher than zero as expected. 
Furthermore, in the = 1/2 case, the curve with no asymmetry seems to be a good fit to 
the experimental data. When both current and resistance asymmetries are considered, 
it leads to a general flattening of the curves, leading to very broad resonances. We have 
chosen k = 0.5 and p = 0.2 as typical asymmetry parameters. We find that increasing p 
to a value higher than 0.2 leads to extremely flattened ib(vd c ) curves and also multiple 
roots in the i>d c (0) curves (next section) which indicates that 0.2 is probably the highest 
resistance asymmetry we can consider in this case. We have also chosen to demonstrate 
trends for k = 0.5 as we feel that for practical SQUIDs this represents one of the highest 
asymmetries. 

5.2. Voltage-Flux characteristics 

Eqn.(17) can also be used to simulate the voltage-flux Vd c — 4> characteristics of a SQUID. 
We have again used Lee et al's SQUID parameter values to simulate Vd c — curves for 
various cases viz. in the absence of any asymmetry in SQUID parameters, in the presence 
of just current asymmetry, in the presence of just resistance asymmetry and also in the 
presence of both current and resistance asymmetries. 

Fig. (4a) shows the va c — curve in the absence of asymmetries. Fig. (4b) shows the 
Vdc — curve in the case when a current asymmetry k = 0.5 is considered. 
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Figure 4. (a) Vd c (4>) curve in the absence of any asymmetry in SQUID parameters, 
(b) i>d c (0) curve in the presence of current asymmetry k = 0.5. The minima exhibit 
resonances and go to zero. In both the figures, bottom curve is for it = 2.4 and top 
curve is for ii, = 3.6. 



Fig. (5a) shows the Vd c — <t> curve in the case when a normal state resistance 
asymmetry p = 0.2 is considered. Fig. (5b) shows the Vd c — <P curve in the case when a 
current asymmetry k = 0.5 and resistance asymmetry p = 0.2 are considered. In each 
of the four cases, normalized bias current 4 ranges from 2.4 to 3.6. 




'0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



(a) (b) 

Figure 5. (a) Vd c {<j)) curve in the presence of resistance asymmetry p = 0.2. (b) i>d c (^) 
curve in the presence of current asymmetry k — 0.5 and resistance asymmetry p = 0.2. 
The minima which indicate resonance positions do not go to zero anymore. In both 
the figures, bottom curve is for % = 2.4 and top curve is for if, = 3.6. 



From fig. (4), we find that in the absence of any asymmetry we have perfectly 
symmetric Vd c {<t>) curves. The difference between the maximum and minimum Vdc 
decreases i.e. Avj c decreases for increasing values of bias current, Avd c goes to zero 
for % = 3.3. Then it again increases for higher bias current values. The complete 
flattening of the Vd c {4>) curve or Avd c going to zero indicates a resonance position. 
With the introduction of current asymmetry (fig. (4)), the voltage values are less than 
the corresponding ones in the absence of asymmetry for a particular bias current; for 
instance, for 4 = 3.6, maximum Vdc in the absence of asymmetry is around 1.62 whereas 
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Figure 6. Comparison between voltage modulation Avdc{ib) where both Avdc and 
% arc in normalized units in cases where no asymmetry is considered and when 
various asymmetries are introduced. The continuous line is the simulated Avdc for 
no asymmetry considered, the dotted line is for the case when current asymmetry is 
considered and k = 0.5, the dashed line is for the case when normal state resistance 
asymmetry is considered and p = 0.2 and the dot-dashed line is for the case when both 
current and resistance asymmetry are considered and k — 0.5 and p — 0.2. 

for the same bias current the maximum voltage decreases to around 1.45 in the presence 
of current asymmetry. Both fig. (5a), where only resistance asymmetry is considered and 
fig. (5b), where both resistance and current asymmetries are considered exhibit skewed 
Vdci.'P) curves. Beyer et al [15] found that t><fc(0) curves got skewed on introduction of 
asymmetry in SQUID parameters. We find that complete flattening of Vdc{<t>) curves is 
absent when resistance asymmetry is introduced i.e. the minima in the curves do not 
go to zero. From Lee's experimental data ( for example see fig.(3a,b)) we find, that at 
% = 3.3, the voltage values for = and = 1/2 are almost the same. This indicates a 
resonance position. Thus we expect flattening of the voltage-flux curves at this position. 
This indicates that even though some current asymmetry may be present in the SQUID 
(see fig. (4b)), resistance asymmetry is probably absent. 

5.3. Voltage modulation AV 

Fig. (6) shows the comparison between voltage modulation Av dc (i fe ) where both Avd c 
and if, are in normalized units in cases where no asymmetry is considered and when 
various asymmetries are introduced. 

We find that the resonance position is almost unaffected when just current asym- 
metry is introduced. However, the magnitude of Av^ decreases by around 20%. On 
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introduction of just resistance asymmetry on the other hand, the magnitude of Ai^ c 
increases, by almost 100% for the first peak and by around 20% for the next couple 
of peaks. The resonance position also seems to shift towards the left (towards smaller 
current values). On introduction of both current and resistance asymmetries, the res- 
onances get rather flattened and there is a general decrease in magnitude of Avd c in 
keeping with the observation for just current asymmetry and the resonance position, 
especially the first peak shifts a little leftward in keeping with the observation for just 
resistance asymmetry. The other peaks do not show any significant shift. 

6. Summary of results 

In this paper, we analytically investigated the effects of transmission line resonances 
on dc SQUID characteristics. We obtain a closed form expression for the ib{vdc) 
characteristics of the SQUID. We use SQUID parameters used by Lee et al [6] and 
compare our theoretical simulations with their experimental data. We find that the 
resonance positions are mainly controlled by IqRs as well as the dielectric constant 
of the substrate €r. We also investigated the effects of introducing asymmetry in 
junction parameters. We find that the resonance positions are not very sensitive to 
the asymmetry. The actual magnitude of the current-voltage curve and AV curves 
are however sensitive to the degree of asymmetry introduced. Our analytical solution 
may provide us with a method of theoretically determining the value of the dielectric 
constant of the substrate if we use that as a fit parameter between experimental and 
theoretical graphs. In an earlier work, we have determined e# of STO thin films using a 
Josephson junction based technique called Josephson Broadband Spectroscopy [16, 17]. 
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Appendix A. 

Enpuku et al's expression for the circulating current J(eqn.(12) of [3]) is given by: 

00 4 

J = Jo(0 2 -9 1 + Y / m n 9 rtn ) + -(f) ex (A.l) 

n=l P 

where m n = M n /L n , M n is the mutual inductance between SQUID and n th resonant 
circuit, L n is the inductance of the n th resonant circuit, 8 r ^ n is the phase corresponding 

to the voltage V r>n across L n , 6 ex is the external applied flux and jo — A I 55 

V \ ~> 2 

n=l 

where (3 is the SQUID inductance parameter and a n is the coupling constant between 
the n th resonant circuit and the lumped SQUID loop given by a n = ^|^rY) • ^ drawback 
in the expression of J is immediately apparent. Substituting the expression for a n , we 
get for J: 

11 00 46 

J = 55 ^ -01 + E m n0 r ,n) + ^ (A.2) 

(1-E«n)^ ^ ^ 

n=l 

There is a problem with this expression in the continuum limit when n — > 00. In this 
limit, 

00 

2 ~ 0l = ~ E m n°r,n (A.3) 
n=l 

=^ J = jj + , which is undefined. 
Appendix B. 

In the following derivation of the SQUID power balance relation, random noise currents 
have been neglected. Multiplying eqn.(12) by 0\ and eqn.(13) by 9 2 , adding the two and 
noting that the voltage across the SQUID averaged over a time T is given by: 

I F 6\dt = v dc = l- [ T 9 2 dt (B.l) 
I Jo I Jo 

we obtain: 

v dc i B = hm 1 / T rft(^ 2 + 61 + 7 (0! - ^ 2 ) 2 - \j{9\ - 9 2 )) (B.2) 

T^oo i Jo 2 
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Eqn.(B.2) can be rewritten as: 

i B v dc = lim 1 f T l-(S 2 + (1 + 2 7 )Z) 2 - j(t)D)dt (B.3) 

T^oo 1 JO Z 

where 5 = Q\ + # 2 and D = 61 — 9 2 . The left hand side of eqn.(B.2) represents the 
input power from the bias current %b- The first two terms on the right hand side of 
eqn.(B.2) represent the power dissipated in the two SQUID shunt resistances, the third 
term represents the power dissipated in the damping resistance while the fourth term 
represents the power dissipated in the resonant circuit. The time averaged voltage as 
a function of bias current and applied flux V(Ib, 3>) for the SQUID can be obtained if 
the voltage waveforms for the two junctions V\ and v 2 are known. When a dc voltage 
Vdc appears across the SQUID, the normalized Josephson currents, sin(6>i) and sin (#2) 
oscillate with time since the relation 



sin(6>j) = sin( J vidt) 



holds, where i = 1,2. Moreover, by averaging over the period of a Josephson oscillation, 
it is simple to show that the fundamental frequency of the Josephson current may be 
given by Vd c /2ii. Therefore, we can express the Josephson current in the presence of Vd c 
as: 

00 

sin(^) = U, n cos{nv dc t + ^ i>n ) % = 1, 2 (B.4) 

71=1 

where %i_ n and ^ n represent the amplitude and phase of the n th Fourier mode 
respectively. The above relation means that the Josephson current can be regarded as 
an ac current generator in the presence of Vdc- The amplitude of the current generator 
is given by [18]: 



I = H = i 2 = pv dc {\ + *4)V2 _ Vdc (B.5) 

The phases of the two current generators are different and the phase difference is given by 
the normalized applied flux as ip\—ip2 = 27T0 . In the present calculation, only the first 
harmonic of the Fourier expansion has been considered i.e. sin(#j) = i^ cos(vd c t + ipi) has 
been used. We have found that this approximation is in good agreement with numerical 
simulations as well over the Josephson frequency range. In this paper, the following 
definition has been used for the Fourier transform: 

X{u) = / X{t)e lut dt (B.6) 
J —00 

and 

X(t) = — / X(uj)e~ tujt duj (B.7) 
2n J-00 

Taking Fourier transform of eqn.(12) and eqn.(13) gives: 
- (1 + x)u 2 PcOi 

= ^i B 6(uj) + l -3 + iu{\ + p% - (1 + ^sin^O + lu-iih - 6 2 ) (B.8) 



On Transmission Line Resonances in High Tq dc SQUIDs 
and 

- (1 - xWPch 

1 l_ _ — _ 

= -i B S(u) - -j + iu(l - p)9 2 - (1 - re)sin(0 2 ) - iui{Oi - 9 2 ) 

Denoting S = 9i + 9 2 and D = 9 X — 9 2) we get: 

a{uS)K + b{u) — {iujp + uj 2 xI3c)D — ?b<5(o>) 



and 

where 

and 

with 

and 



S = 



D = 



{uj 2 /3 c + ioo) 

a(uS) + Kb(uo) - (iuop + u 2 xPc)S ~ j 
(uj 2 (3c + ioj + 2iwy) 

aiuj) = hiu) - I 2 (uj) 
b{u) = h(uj) + I 2 (u) 

iipl 



h(u) = Iir{5{-v dc + iu)e-^ + 5{v dc + u)e 1 ^) , 



I 2 {uo) = l7r(5(-v dc + uo)e~ 1 ^ + S(v dc + u)e^) . 
Using the notation 

s(cu) = (uj 2 p c + iw) , d{oj) = A{uj) — (ou 2 Pc + ioj + 2iwy) , 
r(u) = ipuj + xficu 2 , 

fi = Ka(u) + b(u>) — 5{uj)i b , 

f 2 = P(a{uj) + nb(u)) - Sn<f>A(u)6(u) , 



we get 



(3d{u)h + r{u)f 2 
P[d{u)s{u) - r(u)r(-u)] ' 

-s(u)f 2 + pr^ojfj 
P[d(u)s(u) - r{uj)r{-uj)\ ' 

The inverse Fourier transform leads to 

Ap<t) 



D 



S{t) = 2v dc t - 



-iv dc t 



+ 



2 \d(v dc )s(v dc ) - r(v dc )r(-v dc )} 

+ (l + K)(d(v dc )+r(v dc ))e-^ 



1 - K)(d(v dc ) - r(v dc ))e 
+ c.c. , 



On Transmission Line Resonances in High Tq dc SQUIDs 



15 



I 

+ 77 



-iv dc t 



2 (3[d(v dc )s(v dc ) - r(v dc )r(-v dc )] 

-(lH)(^)-rW)e* 



l-K)((3s(v dc )+r(v dc ))e^ 
+ c.c, (B.21) 



j(t) = const 



-iv dc t 



A{v dc ) 



2 (3[d(v dc )s(v dc ) - r(v dc )r(-v dc )} 

-(l + K)(ps(v dc )-r(-v dc ))e^ 



l-K)(Ps(v dc )+r(v dc ))e^ 



+ c.c. , 



(B.22) 



The constant term in j(t) is given by A(0)(8<f)/[3 — pibf3/4) but plays no role in the 
following. After some tedious algebra it can be shown that (B.3) leads to eqn.(17). 



